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Expressions  arc  derived  for  ttys  notual  cad  self- 
radiatlce  ispedaaces  between  tbe  eleatnts  is  an  array  of 
coexial,  free-flooding,  axially  spaced  ring  transducers 
undergoing  axisyanrbric,  redial  vibratlcs:.  Fluid  clrcu- 
latice  through  the  gape  betirasa  elements  and  around  the 
array  ext  read,  ties  is  accounted  for  in  a  fomilaiioo  of  this 
preblea  which  was  presented  in  a  recent  study  of  single 
*squirters“.  (Cti-foport  An  integral  equation  is 

constructed  which  defines  the  unknown  radial  velocity  dis¬ 
tribution  ct(z)  between  elements  end  on  the  two  sead- 
infinite  cylindrical  surfaces  whit  fa  bracketjfche  array.  It 
is  shown  that  the  radiation  iapelanee  ia  tSe  sun  of  two  cess- 
penents:  (l)  the  irspedaace  evaluated  by  aeans  of  tbe  Math¬ 
ematical  atodcl  crigiaated  by  Robey,  which  assises®  that  the 
active  array  rlaont  is  located  on  an  infinite  rigid  cylin¬ 
drical  baffle,  and  (2)  the  impedance  associated  with  the 
unknown  velocity  distribution  Ot(s),  which  is  therefore  in 
the  nature  of  a  correction  factor  to  Robey’s  impedance. 

Four  aethods  are  presented  for  obtaining  the  function  a(s): 

*  perturbs t loo  solution  which  con  be  used  es  the  sterting 
point  in  a  sore  refined  iteration-type  solution,  a  finite- 
difference  solution  of  the  integral  *  ■'stlon  end  finally  s 
variational  solution.  ^Jtafortunately  *.  variational  solu¬ 
tion  is  less  ettractira^in  tbe  case  of  tbe  array  problem 
than  in  the  case  of  the  single  "squirier”  oscause  only  the 
self -radiation  impedance  Vn  be  obtained  directly  from  the 
stationary  functional  in  terms  of  which  this  approach  is 
formulated.  The  autuel  impedances  wist  be  evaluated  from  the 
uedeternined  coefficients  in  terns  of  which  the  trial  func¬ 
tions  a(s)  are  expressed. 


Dr.  J.  2.  Greenapor.,  of  i  5  Engineering  Heoe&rch 
Associates ,  evaluated  the  inverse  Fourier  transforms  in 
tiaras  of  which  the  solution  ia  prwShted;  an 3  obtained 
quantitative  results  bosh  for  the  single  cylinder  ana  ror 
the  array.  He  vill  present  this  ante rial  in  a  separate 
rtport  entitled  "Axially  Syscstric  Green's  Functions  for 
Cylinders . " 
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(Alternative  subscripts,  vis,  &0>  are  used  to  condense  tvo  equations 

J 

into  one,  tbs  upper  subscript  on  the  left  olds  of  the  equation  being 
associated  vitb  the  upper  signs  and  subscripts  on  the  right  side  of  the 
equation,  and  vice  verse.) 

aean  radius  of  "scuirter”  (Pig.  2) 

a^,So  inner  5cd  outer  radius  of  “aquirter respectively  (rig.  2) 
b  half- width  of  gap  between  eleoents  (rig.  i) 

c  sound  velocity  in  fluid  urdiua 

Gi,'o,Sis,3os,6ia,soa»r,rs,ro,&,As  804  Aa  S^^en’s  functions  and  related 
functions  defined  in  Table  1 

Hq  Hankel  function  of  the  first  kind,  of  order  n.  [This  function  represents 

outgoing  waves  for  the  aasueed  tine  dependence  exp(-lmt)3 

h  half  wall  thickness  of  "squirter"  (Fig.  2) 

J  Bessel  Amctioo  of  order  n 

Q 

k  wave  number.  equtu.  to  a>/c 

i  jj  p  X 

kr  .  idiol  wave  wnsber,  equal  to  (k  -k_)z 

k.  axial  wave  masher 

*» 

I*  half-length  of  array  element  (Pigs.  1  end  2) 

l  half-length  of  array  (Fig.  l) 

p  sound  pressure 

r  jta.'Z  ^yllrnlri— i  cr-JiUijfites 

U  radial  velocity  asplitude  of  active  array  element  (Figs.  3  and  k) 

u(z),u.(z),u  (c)  radial  velocity  distribution  over  cylindrical  surfaces  (r=a), 
(r-^)  and  (r=aQ)  respectively,  positive  outward  (Figs.  2  and  3) 


uB(s),uo(s)  syET^atric  (even)  and  antisymmetric  (odd)  component  of 

°  velocity  distribution  over  cylindrical  surface  (r«=a)  (Fig.  k) 

Z..  self-radiation  lapedsnce  of  crray  element  j,  in  units  of  forcc/velocity, 

e<3ual  *>  r  *  <Zj*U 


Vr 

V* 


mutual  radiation  impedance  experienced  by  array  element  k  as  a  result 
of  vibration  of  array  element  J ,  In  units  of  force/velocitv  “'"cal  to 

«V  *  * 

radiation  impedance  obtained  from  Hooey 1  s  mathematical  model  for  which 
a(s)»0  (transducer  elements  located  on  infinite,  rigid,  cylindrical  baffle) 

correction  factor  associoted  with  velocity  distribution  a(z),  to  be 
added  to  (2jfe)r 


*Otbsr  symbols  arc  defined  in  the  text. 


sf  region  of  z-axis  consisting  of  the  grp«  between  am*  elsawssa ,  sad 
of  the  two  seai-iilfinite  rrlitirical  surfaces  (z  >  l,  %  <  t)  vMch 
brfe&et  the  array  (fig.  3} 

sf+  positive  half  of  z? 

tg  location  of  aldplanc  of  gaps  between  neighboring  array  elements  (fig.  l) 

Zj  locetion  of  aidpisae  of  active  array  element  (fig.  4) 

n.  location  of  aldplane  of  array  element  (fig.  l);  la  iaqpeisnce  calcula¬ 

tions.  location  of  Inactive  element  (fig.  4) 

a(t)  radial  velocity  in  the  region  a.  (fig.  3)  normalized  to  velocily 
amplitude  V  of  active  element  •  ar  array 

a :  (z),a (z)  symmetric  (even)  and  antisymmetric  (odd)  component  of  velocity 
distribution  in  region  z,,  normalized  to  U/2 

8.  ,6  coefficients  In  the  relation  of  the  velocity  0  of  the  naan  transducer 
surface  with,  respectively,  u  and  u  ,  defined  in  Bq.  H.3,  *  1  for 
snail  h/a  »  1  0 

fi(z-£/)  Dirac  delta  function,  f a(s/)  dz*  =  Ct( 2) 


s  Poisson's  ratio  of  transducer  arterial 

p  density  of  fluid  medium 

}(r,z)  velocity  potential  (outward  velocity  is  -bt/br);  subscripts  "1“  and 
"o"  refer,  respectively,  to  regions  r  ^Caj^end  r  >  ao;  subscripts  "0" 
and  'a*  refer,  respectively,  to  symmetric  (eve n)""aad  antisymmetric 
(odd)  components  with  regard  to  z 

m  circular  frequency  [with  this  notation,  s  massive  reactance  is  negative; 

harmonic  tine  dependence  factor  exp(-iajt)  which  multiplies  the 
velocities  and  the  potentials,  has  been  suppressed  throughout  this 
report] 
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I.  Scope  and  Background  of  This  Study 

The  purpose  of  this  study  is  to  evaluate  the  radiation  impedances  for  the  el¬ 
ements  of  an  array  consisting  of  identical  free-flooding  coaxial  ring  transducers 
vibrating  in  their  radial,  axis  ymrne  trie  mode  (Figures  1  and  2).  The  transducer 
elements  are  separated  by  equal  axial  distances,  thus  permitting  the  acoustic  medium 
to  flow  both  axially  at  the  ends  of  the  array  and  radially  between  elements.  The 
lapedances  take  into  account  the  amfeined  radiation  loading  acting  on  the  inner  and 
outer  faces  of  individual  transducer  elements,  when  only  one  element  in  the  array 
is  active.  The  result  of  this  study  will  thus  be  in  the  fora  of  self-  and  mutual 
lapedances  of  the  array  eleesnts.  Cogitations  will  be  performed  in  the  resonant 
frequency  range  of  the  transducer  elements,  the  variable  parameter  being  element 
spacing.  Even  though  this  study  was  stimulated  by  the  promising  experimental  re¬ 
sults  recently  obtained  by  the  U.S.Kavui  Research  Laboratory  on  an  array  of  mffx- 
tostrictlve  ring  transducers,  the  present  aiwlysis  la  not  limited  to  magnetostric- 
tive  transducers  nor,  in  fact,  to  free-flooding,  axially  spaced  elements.  Wien 
conbined  with  the  analysis  of  Individual  free-flooding  or  solid  cylindrical  radia¬ 
tors  presented  in  an  earlier  report,1  the  present  analysis  can  be  applied  to  the 
evaluation  of  mutual  lapedances  of  ring  and  piston  elements  on  s  cylindrical 
baffle  of  finite  length. 

The  array  analysis  can  be  followed  independently  of  reference  1  if  the  reader 
is  willing  to  accept  without  proof  certain  results  derived  therein.  It  was  shown 
in  this  earlier  study  that  the  radiation  iffpetlances  of  an  isolated  transducer  ele¬ 
ment  can  be  expressed  as  the  sua  of  two  components: 

Z  »  Zr  ♦  ZQ  (T.l) 

2  14  4 

is  the  impedance  computed  by  Robey,  Greenspon,  and  Sbercon  from  Robey's 
cat  heretical  model  which  is  frequently  used  for  cylindrical  transducers.  This 
model  assures  that  the  radiating  surface  is  extended  to  infinity  by  rigid  cylin¬ 
drical  baffles  which  prevent  circulation  of  the  fluid  around  the  transducer  edges.9 
This  assumption,  in  combination  with  a  Green's  function  constructed  so  as  to  have 
its  radial  derivative  cG/ar'  vanish  on  the  outer  transducer  surface  (r'«a). 


•The  radiation  loading  cm  the  inside  of  a  free- flooding  cylindrical  transducer 
was  also  studied  by  Robey  by  Beans  of  two  different  code la ,  which  unfortunately 
also  place  restrictions  on  the  fluid  flew  around  the  transducer  edges.8’7 
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eliminates  the  seed  of  solving  an  integral  equation  to  obtain  the  potentials. 

With  this  modal  an  expression  for  the  far  field  potential  55a  ba  obtained  ana¬ 
lytically  in  closed  fora.'*  Jbe  impedance  ZQ  in  S3.  1.1  is  in  the  nature  of  a  cor¬ 
rection  factor  associated  with  the  radial  flow  Telocity  distribution  a(z)  across 
the  cylindrical  surface  extending  the  transducers.  is  in  the  fora  of  an  inte¬ 
gral  vboee  integrand  is  proportional  to  a{s).  The  function  a( z)  satisfies  u  non- 
hoBogeneous  Fredholn  integral  equation.  Unfortunately  a  direct  solution  of  this 
equation  is  not  practicable.  The  analytical  effort  is  therefore  neatly  directed 
at  approx iaating  the  unknown  function  a(z). 

In  ref.  1  a  variational  technique  t*s  developed  which  parallels  the  Lerine- 
Scbwlnger  varistionol  procedure  for  confuting  scattering  cross  sections  in  dif¬ 
fraction  problem.*  In  this  procedure,  a  functional  Jja}  ia  constructed  froa 
the  integral  equation.  This  functional  can  be  shewn  tc  be  stationary  with  respect 
to  first  or'ler  variations  about  the  solution  a(z)  c'"  the  Integral  equation.  If 
rctio  of  transducer  wall  thickness  to  radius  Is  swell,  the  stationary  value 
of  J[i]  is  proportional  to  iQ.  A  Rayleigh-Bitz-type  oilcaUticn  is  used  to  cos- 
pute  Z  ,  starting  with  a  trial  function  for  a(z).  The  met  cobboo  use  of  the 
Rsyleigh-Ritz  oethod  is  the  evaluation  of  the  natural  frequencies  of  a  vibrating 
systen.  These  frequencies  are  toe  roots  of  the  determinant  of  the  coefficients 
in  o  set  of  snaulteneous  homogeneous  equations  derived  frees  a  variational  prin¬ 
ciple.  The  retural  frequencies  cesi  therefore  be  evaluated  without  having  to 
compute  the  undetemlned  coefficients  in  the  sssuaed  aodai  configuration  of  the 
vibrating  body.  In  the  saaljolB  of  the  single,  thin-walled  "squlrter,"  ZQ  can  be 
siallorly  obtained  without  jeeviously  solving  for  the  unknown  coefficients  In 
terns  of  which  a(z)  bss  been  expressed.  Unfortunately,  only  its  self-radiation 
ispedaace  can  be  deterained  in  this  fashion  when  the  transducer  is  on  elesent  in 
an  array.  In  order  to  evaluate  the  nut  uni  izpedances,  the  eet  of  sisailtaneeus  al¬ 
gebraic  equations  derived  froa  the  variational  principle  oust  be  solved  for  the 
unknown  coefficients  in  the  trial  function  o(s).  Thus,  while  in  the  single  trans¬ 
ducer  prcblen  the  variation  a  ssothod  is  considerably  cove  efficient  ttss  =Ucrz»- 
tive  methods  for  solving  «».•  integral  eqo-iion,  this  io  not  the  =:<-»  array 

problea. 

Ve  will  therefore  present  several  other  techniques  for  dealing  with  the 
Integral  equation.  Eerier  ce  with  mzaarical  calculations  will  show  which 

~~  Vl  bibliography  of  variational  analyres  of  diffraction  problem 

is  given  in  the  bibliography  of  ref.  1. 


procedure  It  nost  efficient  for .»  glen  conbiustloa  of  array  paraatters.  One  of 
the  four  technique*  bere-  preecntei  involves  finite-difference  calculations.  Evan 

'*\  /=  V  jt  N*  "8'0 

though  such  calculations  here  recently  been  usad  by  various  .authors  ’  to  solve 
the  Selafcoltz  integral  equation  foraulated  la  terns  of  -tt*  frae-space  Green’s 
function,  their  forsalatian  of  the  prbblsn,  is  basioally  different  fro*  the  present 
one.  In  references  8  sad  9,  the  unknown  function  in  the  Integral  equation  is  the 
potential  at  the  solid  surfaces-  vhitii  constitute  the  boundaries  erf  the  acoustic 
aediuK.  Za  the  esse  of  an  arrsy,  these  surfaces  would  include  all  inactive  as 
well  as  active  array  elenents.  In  the  analysis  presented  here,  the  unknown  func¬ 
tion  is  a(r),  the  radial  velocity  distribution  in  the  annular  specs  between  the 
ring  elements  and  in  the  two  seal-infinite  anoulrar  specie  extending  the  array, 

Thin  unknown  fun?  »  a(z)  ashes  a  less  important  contribution  to  the  radiation 
ispedance  then  does  the  unknown  potential  over  the  transducer  surfaces  in  the 
above -auctioned  analyses.  A  cowparvble  accuracy  in  the  iapedance  calculations 
can  therefore  be  expected  froa  the  present  approach  using  s  coarser  finite- 
difference  spacing. 

Quantitative  results  will  he  presented  in  a  separate  report  by  Sr.  J.  E. 
Greenspon,  of  J  0  Engineering  Besaarch  Associates,  entitled  "Axially  Synaetrie 
Green’s  Functions  far  Cylinders."  C*.  Greenspon  Is  obtaining  numerical  results 
for  the  transducer  peraneters  corresponding  to  the  ICS.  array. 

H.  Foraailatlon  of  the  Problea 

Following  Bobey’o  exaaple  the  potentials  are  expressed  in  terse  of  two 
Green's  functions  whose  radial  derivatives  vanish  on  one  or  the  other  of  the  two 
infinite  concentric  cylindrical  surfaces  on  which  the  inner  and  outer  transducer 
faces  lie.  Two  potentials  are  thus  constructed  *0,  in  the  outer  region  r  >  a^, 
and  I,  in  the  inner,  cylindrical  region  r  <  a,^  (Fig.  2).  These  potentials  are  ex¬ 
pressed  in  terse  of  the  radial  velocity  distributions  u^z)  and  uq(z)  over  these 
two  coaxial  cylindrical  boundaries:* 

*5 

4^(r,s)  »  t  2sa^  j  u^(z')  &^{r.a^,  z~z')  d z'  (H.l) 

"Sere  and  elsewhere  in  this  report  alternative  subscripts  and  signs  have  been 
used  to  condense  two  equations  into  one,  the  upper  subscript  on  the  left  side  of 
the  equation  being  associated  with  the  upper  signs  and  subscripts  on  the  right 
side  of  the  equation,  and  vice  versa. 


B*  velocity  ittlibM  over  Inactive  transducer  nrftci  cod  equals  a  known 
constant,  »y  U  in  -the  regiai-ebere  the  active  transducer  la  located.*  the  ana¬ 
lytical  difficulty  lias  1b  the  Sect  that  the  radial  velocities  are  known  only  on 
the  transducer  surfaces  (H|.  3).  Before  defining  the  Green's  function*  it  is 
cocrenient  to  express  the.  velocities  u^Cs}  *nd  uq(s)  in  term  of  the  velocity 
u{z)  of  the  men  surface  (r*a).  for  this  purpose  ~a  asiua*  an  jncoapreasible 
potential  in  the  annular  space  <  r  <  aQ.  This  laplies  that  the  ratio  pf  mil  - 
thickness  to  acoustic  wavelength  is  negligible.**  With  this  assumption,  the  bal¬ 
ance  of  mss- inflow  anu  outflow  yields  the  following  relations  between  velocities: 


u^(s)  *  u(z)(lt  ~)“A,  for  z  in  region  zf. 


£vl 


(n,2s) 


the  region  zf  encompasses  those  portions  of  the  annular  space"  (a^  <  r  <  aQ) 
which  contain  acoustic  fluid,  rather  than  transducer  elamnts.  She  region  zf  thus 
includes  the  gaps  between  the  array  elements  end  the  two  seal-infinite  annuli 
s  <  £  and  z  >  t,  which  extend  the  array: 


-l  — ,  zg+b  • 

%  a*  g  *g-t>  * 


dz 


where  s_  is  the  coordinate  of  the  plane  of  symetzy  of  a  gap  between  two  nelgh- 
e 

boring  array  eleaeuts 
coordinates  of  the  plane  of  sy 
Sable.  The  velocity  of  the  inner  and  outer  transducer  faces  are  related  to  the 
velocity  ef  the  seen  transducer  surfs cs  ns  follows: 


She  values  of  z  are  given  in  the  Sable  in  Fig.  1.  The 
try  of  the  elements  are  also  given  in  this 


«  (z)  -  u(z)B  (1+  I)'1 

?  ? 


for  Zj-L  <  z  <  Zj+L 


<n.2fb) 


where  the  coefficients  9  esbody  the  Poisson's  ratio  of  the  transducer  wterlal: 

for  piezoelectric  transducers 


s„  ■  (u  |)(»  r> 


for  aagnetostrlctive  transducers 


Cn.3) 


"She  tim  depandenee  exp(-iajt)  has  been  suppressed  throuf^aout  this  report  for 
the  sake  of  brevity. 

**£t  was  pointed  out  in  raf.  1  that  compressibility  of  the  fluid  in  the  annular 
region  could  be  included  is  the  analyrls  but  that  this  would  result  in  tw:  coupled 

IsUay*!  cyretiznz  isztusd  vt  a  aizigle  ifivegril  awuSaxuci.  lu  the  case  sf  tiaa  tfhli 

masatoitrictirs  array,  the  ratio  of  well  thickness  to  acoustic  wavelength  is  0.07. 
She  aacusptice  of  an  incompressible  potential  is  therefore  Justified  within  the 
licit £  of  eccuracy  required  for  design  purposes. 
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The  subscript  4  Kill  be  used  for  the  acfctra  transduca?,  sad  the  tube  dipt  k  for 
the  las ct ire  cransdactrs.  the  origin  of  the  x-coordiasrfce  is  located  Is  the  plane 
of  ayitry  of  the  array. 

She  redial  Telocity  distribution  over  the  seen  cylindrical  surface,  (r^s),  ia 
expressed  as  follove  (see  fig.  3}: 

u(s)  *0  on  actlre  array  eleeeot;  *j-Ii  <  *  <  s^+L 

»  U  a{z;  for  z  is  region  z, 

=  0  os  inactive  array  elenests,  z^-L  <  z  <  z^+L  (see  Table 

in  fig.  1)  (U.k) 

The  unknown  function  a(z)  Is,  In  general,  ccaplex. 

The  Green's  functions  for  the  two  Integrals  In  £9.  U.l  are  In  the  for*  of 
inverse  Fourier  transfer**: 


00 


in 


where  the  transforms  cf  the  Green's  function  are 

H  (k  r) 

Go(r»°ojks)  "  aS"  SrH1(krao) 


(n.5) 

(11.60) 

(H.6b) 


It  can  be  verified  that 

<JG_/dr'  J  «  0 

n*  ^ 

i 


v^. 


Tbs  potent  is  la  can  be  obtained  by  evaluating  the  inverse  Fourier  transferns, 

Eqs.  II. 5,  and  then  performing  the  s'- integration  Indicated  in  Bq.  U.l.  Alter¬ 
natively,  the  order  of  lntegraticc  can  be  reversed.  When  Bqs.  II. 5  are  (substi¬ 
tuted  is  5q.  II. 1,  the  z'-integratioa  can  be  performed  fr?st,by  asking  use  of 
the  definition  of  the  Fourier  transfer®  of  the  velocity  distribution: 


9 

»(k_)  *  f  a(z')  cxp(-iX„i‘” )  dz# 
“  J 


(n.8) 


5 


/ 


She  potentials  in  Bq.  H.l  can  then  he  written  in  the  foam  of  «n  inverse  Fourier 


transform 


*oi->TJ  ac  f  c(iE)*xp{ikEz) 


(2X.S) 


As  tiie  dsvelopnent  of  this  analysis  need  not  he  specialized  to  O’M  or  the  other 
of  these  procedures ,  their  relative  convenience  can  he  coapared  after  acre  egc- 
twsrieace  has  been  obtained  in  the  evnlnatica  of  numerical  results. 

Wh«  ve  substitute  Bqs.  U.2  and  k  in  Sq.  H.l.  we  can  express  the  potentials 
generated  by  tur  vibrations  of  the  Jth  eleacnt  of  the  array  as  follows: 

#  (r,t)  *t2xa  0iBo  f  3„(r,a  ,z -z')dz'  +  I  a(z')G0(r,eo,z-a')dz']  (H.10) 
i  1  z,-L  11  z  1  1 

st  * 

The  unknown  velocity  distribution  a(z')  in  these  integrals  is  defined  by  the  re- 
quireaent  that  at  every  location  the  pressure  differential  between  the  concentric 
cylindrical  surfaces  (r=z, )  end  (r=«o)  balance  the  radial  cceponent  of  the  iner¬ 
tia  force  exerted  by  the  fluid  in  the  annular  region  between  these  tw>  surfaces. 

The  radial  acceleration  averaged  over  the  radial  thickness  of  the  annulus  is, 
taken  positive  outward 

u(a)  =>  -  irOa(z)  for  s  in  region  zf  (H-ll) 

The  racial  coapaaent  of  the  inertia  force  associated  with  a  differential  voluae 
decent  Sba&pdz  of  the  fluid  annulus  Is  therefore 

4P(z)  =  2IaxJa(z)p  hud?  dz  (H.12) 

The  outward  force  associated  with  tha  pressure  differential  is 

d P{z)  «  [SjP^tr.)  -  a^Cz)]  ds?  dz  {H.13a) 

In  turns  of  the  potentials  and  of  the  oean  radius  sod  shell  thickness  this 
baccaas 

4F(t)  =  -fcspa{(l-  -  (l+  ®  ( II.  13b) 

Whan  we  set  the  sun  of  Eqs.  n.12  and  13b,  equal  to  zero,  as  required  by 
d’Alezd>ert*o  principle ,  we  obtain  the  following  relation  between  the  potentials 
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and  the  unknown  Telocity  distribution  a{z): 

(1+  f)*0(®Q>2)  -  (l-  f)*1(«1,s)  3  '2h0a(s>  for  z  in  region  zf  (U.14) 

For  tvo  of  the  configurations  analyzed  in  reference  1,  viz.  the  infinitely  thin- 
vslled,  cylindrical  radiator,  and  the  solid  cylinder,  the  right  hand  side  or  this 
equation  vanishes  end  Eq.  H.lb  reduces  to  a  continuity  requirement  of  and  *0 
on  the  surface  (r=a).  When  ve  substitute  the  expression  far  the  potentials, 

Eqa.  H.1C-,  in  Eq.  U.lU  we  obtain  the  non-hoaogeneoua  Freihola  equation  which 
defines  unknown  function  a(z): 

z^+1 

j  [r(s-z/)+  ~6{z-z/)}a{z,)dz#  »  -  f  d(s-5#)d* * ,  for  z  in  region  zf  (n.15) 

sf  v1. 

In  this  equation  the  oyaocls  T  and  A  have  been  used  to  represent  two  linear  <xm- 

Mirations  of  the  Green's  functions  Q,  and  G  : 

J.  c 

r(z-z')  -  (1+  |)  Oo(ao,ao,z-z')  ♦  (l-  |)  G^a^.z-z')  (n.lfe) 

b{z~z'}  -  8o(l+  |)  Go(ao,ao,z-z')  ♦  #£-  f)  G^a^.z-s')  (H.lfb) 

Having  thus  derived  the  integral  equations  vhich  a{z)  must  satisfy,  we  proceed 
to  express  the  iapedances  in  terns  of  this  sane  velocity  distribution  a(z). 

HI.  Badjatioa  Inpedance 

The  radiation  loading  exerted  by  the  Jth  elesaent  of  the  array  on  the  kth 
elsssent  is  obtained  by  integrating  over  the  surface  of  the  latter  the  pressure 
differential,  Eq.  II. 13b,  where  the  potentials  are  associated  exclusively  with 
the  notion  of  the  Jth  transducer  eleaoat: 

V  ■  -  ^  J*  Id  *  |«0(-0-‘)  •  d-  i'lqf.,.--)!  *  <m.i) 

VL 

A  radially  inward  directed  radiation  load  is  positive  in  this  notation.  When  ve 
substitute  the  expression  for  the  potentials,  Eq.  H-10,  and  coking  use  of  the 
abbreviations  for  cot&inationa  of  Green's  functions,  Sqs.  H.16,  ve  finally  obtain 
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the  following  expression  for  the  radiation  i^edanca 
*j+L  rS*^ 

x  .  {2xa)2Jfflpf  tf  A(*-s#)  da*  +  fa  (s')  r(sV)  da']  ds  (XXX.2) 

VL  sf 

As  in  the  case  of  the  self-impedance  of  the  single  transducer,  the  sutual  iaped- 
snese  ors  5vis  to  be  the  sux  of  two  components:  (l)  the  impedance  obtained  frca 
Bobey's  aatheoetlcal  aodel: 

a^+L  Sj+L 

(2Jk)r  -  -(a»)2ifflpf  f  d(s-z')  ds'  ds  (m.3) 

V*  VL 

and  (2)  an  ispe dance  ccapooent  associated  with  the  uckaoim  velocity  distribution 
a(s): 

Zj+1 

-  -(2*a)2i*p  f  Ja(x') r(s-s')  ds'  ds  (HI.*) 

VL  2f 

The  problee  is  thus  to  find  the  solution  c(s)  of  the  integral  equation,  Bq.  11.15 
and  to  substitute  this  solution  in  the  above  expressions  far  the  radiation 
impedances .  The  velocities  of  the  individual  transducer  elements  can  then  be 
obtained  frcn  the  usual  array  analysis.  Finally,  the  velocity  distribution  a(z) 
can  be  used  to  obtain  the  far  field  potentials  fro*  the  expressions  given  for  the 
single  "squirter"  in  Sections  VII  and  Vlll  of  ref.  1. 

IV.  Fonsulatlon  of  the  Problea  In  Terns  of  Sy— etric  (Bren)  and  Aatiswetric 
fold)  Velocity  Distributions  “““ 

la  Robey's  saatbeaBtlcal  aodel  tbs  velocity  distribution,  and  hence  the  poten¬ 
tial,  are  always  sysnetric  about  the  plane  of  syanetry  of  the  active  transducer 
ring.  In  the  present  aodel,  which  allows  for  fluid  flow  between  elements,  this 
is  not  the  case  except  when  an  active  elesssnt  is  located  at  the  center  of  the  array. 
In  two  of  the  four  solutions  presented  in  Sections  V  to  VII,  it  is  advantageous  to 
consider  the  velocity  0  of  the  active  array  elenent  e~  the  resultant  of  two  cocpo- 
nent  velocity  distributions ,  one  of  which  is  tsyrse-  ~  !-•  about  the  plane  of  synwtry 
of  the  array,  (z=0),  sad  the  other  cntisyaaetric.  «s  seen  froa  Fig.  b,  this  anounto 
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i 

I  to  considering  the  Telocity  0  of  the  active  element. 


u(z)  a  U  for  Sj  -I  <  2  <  Sj  +  Ii  (17.1a) 

as  The  superposition  of  (l)  the  Telocity  distribution  associated  vith  two  ele- 
aents  located,  respectively,  at  z^  and  -z^,  excited  in  phase,  with  half  the 
original  amplitude 

us(s)  =>|  for  +Sj-L<i  <iz^  +  Xi  (iV.lb) 

and  (2)  the  Telocity  distribution  of  the  ease  two  elements  excited  oat  of  phase 
ua(z)  =  ^  for  Zj  -  L  <  z  <  a ^  +  !• 

ua(z)  =•  -  ~  for  -Zj  -  1  <  2  <  -Zj  +  L  (lY.lc) 

For  a  centrally  located  element,  the  antisymmetric  velocity  distribution  does  not 
arise,  and 

u(z)  *  us(z)  =  0  for  -L  <  z  <  1  (zy*0)  (iv.ld) 

Because  of  the  linearity  of  the  problem,  we  can  solve  independently  the  symmetric 
case.  Eg.  IV.  lb,  which  gives  rise  to  a  sysvetric  potential. 


(iV.ga) 


and  the  antis yroetric  case,  Bq.  JV.le,  which  gives  rise  to  an  antisymmetric  po¬ 
tential, 

4n(r,z)  *  -  5a(r,-z)  (I7.2b) 

The  unknown  velocity  distributions  corresponding  to  these  two  potentials 
batiafy  ainJlar  syssatry  relations: 


•  3  (-S> 

0  O 

«c(-)  -  -«,(-=)  (IV. 3) 


It  is  advantageous  to  aorcalize  the  two  coop-nents  of  the  unknown  velocity  dis¬ 
tribution  to  D/2.  With  this  convention,  usd  using  Bqa.  IV. 3,  tie  velocity  dis¬ 
tribution  la  region  z,  can  be  expressed  as 

u(z)  »  |  (as(z)  +  o8(*)J  for  z  >  0 

>  in  region  zf 

«(a)  •  |  lae(s)  -  aa(z)3  for  a  <  0  j  (rv.b) 

/ 
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It  is  6mb  frca  I«a.  IV.£  sal  3  that  if  i  cod  a  satisfy  the  integral  equations, 
ft}*.  11.14  25,  in  tbs  positive  half  of  the  region  i<.  they  automatically  satisfy 

these  equations  is  the  negative  ho  If  of  sf.  She  sum  applies  to  and  aft.  Vt  «m 
therefore  confine  the  integral  equation  to  the  poeitive  belf  of  the  region  zr,  which 
will  he  designated  henceforth  by  the  symbol  sf+.  The  advantage  in  splitting  the  po- 
tentiel  Into  sy*ertriq  end  sntiayneetric  ccapooeut*  results  free  the  feet  that,  with 
the  finite-difference  and  variational  techniques  described  below,  it  la  less  labo¬ 
rious  to  satisfy  two  separate  integral  equations  in  the  subregion  zf+,  than  to  sat¬ 
isfy  coe  integral  equation  ore?  all  of  the  region  -+• 

Ifce  resultant  potential  generated  by  the  active  trandueer  located  at  3a 

#(r,s)  *  *a{r,z)  ♦  *a(r,z)  (lV»5a) 

which,  in  view  of  ft} s.  IV. 2a,  can  be  written  aa 

•=  »g(r,+  |z|)  t  ^(r/flzt)  (lV.5b) 


Ve  need  therefore  only  evaluate  the  two  component  potentials  for  positive  values 
of  z  to  obtain  inforaation  on  the  resultant  potential  over  the  whole  rengs  of  z. 

Restriction  of  tto  analysis  to  poeitive  values  of  s  does  net  Intar fere  with 
obtaining  all  of  the  desired  radiation  impedances .  S»  actual  impedance  is  deter- 
ained  only  by  tav  separation  |zj-zk|  between  the  aetivs  and  inactive  array  element. 
The  self-ispedanct  is  the  k\s*e  for  eieaaats  sjnawstricelly  located  at  Zj  end  -Zj . 

Ve  can  therefore  rest.-*  t.tt  analysis  to  active  elessents  for  which  >  0,  but  z^ 
can  cf  course  be  negative  as  well  as  positive. 

wc  atoll  row  derive  the  co*5«ncat  Green's  functions  which  are  applicable  to 
the  sysestrtc  and  antisysaetrlc  potentials.  The  Green's  functions  in  Bq.  0.5  can 
be  expressed  in  teres  of  an  inter  »e  Fourier  cosine  transform  by  noting  that  the 
imaginary  component  of  the  expc&sntial  doe3  not  contribute  t&  the  value  of  the  in¬ 
tegral,  because  the  Green's  function  is  even  in  Xz: 


G9(r»Vs-s'> 
I  i 


cos[k.(s-z')J  <ft 

5  S 


(IV.S) 


Expanding  tt»  coains  in  the  integrand,  the  Green's  function  cau  be  written  as  the 
sues  of  a  cosqsccent  which  is  even,  or  sysssctric  about  both  z=0  and  z**0,  and  of  a 
coe^KSwat  which  is  odd,  or  ontieysaetric: 


So(r,s0j6.£'} 

i  i  8 


(XV.7c) 


f 

sjj,*v,z-z\  «  ij  5o(r,«j,ks?  •inCkE*}«i£&s*/i  dk. 


(BT.Tb) 


If  both  the  Green's  function  ana  the  velocity  distributie  ,  la  the  integrand  of 
Eq.  II.  1,  are  express a-i  as  a  bus  of  epaastric  and  antiayaastric  teres,  the  cross 
terss  are  found  cot  to  contribute  to  the  value  of  the  integral: 

CO 

±2**j  J  ty*')s  8|{*',«j..s-=#)a+^(ii#)Ei^(r1aJs-s/)a3an/  (IV. 8) 

<%(• 

The  two  products  in  the  integrand  thus  correspond  respectively  to  #  and  i  . 

&  ft 

Since  the  contribution  of  the  integral  to  these  two  notentiale  is  the  saae  for 
the  positive  and  negative  halves  of  the  c'-axic,  these  two  coapoaents  can  finally 


be  written  as 


00 

yr,*)a  “  -  «wr?  J  S^r.a^.a-aOg 

C* 

Jo(r,-)Q  -  t  4«a0  j  UoCs'Jg.  SoCr*a0*z*2,)8  4-*' 
i  i  ^  i  li 


(IV.9») 


(iv.jb) 


Hie  integration  over  z'  can  be  performed  before  the  k^- integration  indicated  in 
Eqn.  IV. 7  by  making  uar  of  the  Fourier  transform  of  ug  and  ,  respectively: 


CO 

*o0->*)8  “  i  2aJVr'° o’V  u{ks)s  008  V  **! 


(IV. 10») 


tov*»'/8  -  -  —  J'  ’  J*"*'  —  "s"  t  - - 

0 

<B 

^(r;£/B  *  -  2  a-  f°j(r»*o>,t*)  u(ks)8  ks»  <&s  (XV.lCb) 

*  6  * 

She  linear  ccBhinatiens  of  Green's  functions,  Eqs.  11.16  used  in  the  integral 
equation.  Eq.  11:15  c£n  sou  also  b=  split  into  syaaetric  and  antisymnsstric  coa- 
ponents: 

r|(=_E  )  *  (1+  a^5o^°o’9o,s’s,^g  +  e^  6i^ai,,ll'*“a#^g  (IV. lie) 


-  eQ(i+  +  SjCi-  (w.ub) 

She  two  ,&&q oeeat  velocity  distributions,  Or  end  a  satisfy  tbs  two  uncoupled  in- 
tegrsl  fttjcfls tiona  which  take  the  piece  of  >9.  H.15: 

f 

J  [r  (a-s#)  +  ~^(z«s#)ta  (a,)dz/  *  -  J  htz-z*)  ds*  for  a  in  region  zf+ 

Sa.  8  •  •  f  * 


3  (l?.is>) 

Robey’s  radiation  inpeftance,  Bqs.  HI. 3,  now  becomes 

z^+L  zytL 

(ZJk)r  -  -(2»)2iop  J  J  ti8(z-s')  +  da(s-z4)}  dr'ds  (IT.13) 
VL  V* 

The  correction,  factor  associated  with  the  radial  velocity  distribution  a(s)  is 


(Zjk)a  -  -(Sxafiapj  J  (osa(z4)rs{z-z')  ♦  aa(s')ra(*-t')3  dt'dz  (iV.lk) 

VL  •» 

If  cither  or  z^  are  zero  than  the  antisymmetric  teraa  and  d#  drop  out  of 
the  integrands  la  Bqs.  IV.13  and  lb.  If  all  the  eleaects  in  the  array  are  identi¬ 
cal  and  driven  with  the  stae  signal,  the  velocity  distribution  and  banes  the  far 
field  potential  do  not  contain  any  antisymmetric  components.  We  will  now  diecuso 
four  approaches  to  the  solution  of  the  Integral  equation. 


V.  Perturbation  and  Iteration  Solution  of  the  la 


la  these  two  epproches  nothing  is  gained  by  dealing  separately  with  syc- 
aetric  and  antisymmetric  velocity  cragwoenta.  The  crudest  of  the  four  approaches 
to  be  described  is  the  perturbation  solution  in  which  we  actually  circurerent  the 
need  of  solving  tbs  integral  equation.  Ibis  approach  starts  with  the  near  field 
potentials  associated  with  Bobey*#  aodal.  They  are  computed  fro*  Eq.  II. iu  where 
a(s')  is  set  equal  to  zero.  Ceslsasting  tbsis  unperturbed  potentials  by  and 
substituting  then  in  Sq.  I2.lt,  we  can  solve  directly  for  the  perturbation  solution 
for  a(s) 


When  ue  substitute  tha  uaperturbad  potential*  w  computed  from  *5*.  H.  10,  where 
a{z)  tes  been  set  equal  to  *ero,  «ad  using  tte  Qrcau:*  function  combination  la 
Sq.  H.lSb,  Sq.  Y.1  can  then  he  vrlvftee  »re  explicitly  u: 


o(0)(s) 


d(z-z')  da# 


(V.8) 


Xu  the  component  subregion*  of  zf  which  correspond  to  the  gape  between  the  array 
elements,  a  somewhat  more  refined  solution  for  a(z)  can  be  obtained  if  one  defines 
an  effective  radial  annular  width  equal  to  the  transducer  well  thickness  (Sh)  in* 
ertnsed  by  u  length  £h  embodying  the  accession  to  Inertia  of  a  silt  In  a  baffle. 

The  Introduction  of  an  accession  to  Inertia  la  not  useful  in  any  of  the  other  ap¬ 
proaches,  where  the  potentials  tQ  and  themselves  eabody  this  reactive  import¬ 
ance  CGOpOSliUt. 

Since,  for  design  purposes,  we  probably  need  a  more  accurate  expression  than 
the  perturbation  solution,  we  use  the  value  Rj.  V.2,  as  a  start  for  an 

iteration  procedure.  The  first  step  consists  in  computing  values  for  by 
substituting  the  perturbation  solution  in  Bq.  H.10.  When  we  substitute  the 
potentials  thus  computed  in  the  integral  equation,  5q.  11.15,  we  obtain  the  follow¬ 
ing  score  refined  expression  ?cr  a(z): 

/,+L 

d(z-z')  dz'  ♦  J  r(z-z')o{0V)  dz')  (V.3) 

Zfl  sf 

This  iteration  can  be  repeated  until  the  fluctuations  of  successive  solutions 
a(s)  are  deemed  sufficiently  small.  Convergence  of  this  process  can  best  be  test¬ 
ed  eapirically  by  psrforuing  the  actual  numerical  calculations  and  comparing  the 
results  of  seesosaive  Iteration  with  results  obtained  from  the  other  approaches  to 
the  solution  of  the  integral  equation. 


VI.  Pialte-Plfference  Solution  of  the  Integral  Sanation 

In  this  approach,  the  region  t,  is  divided  into  a  ouefcer  Q  of  subregions. 

Each  subregion  is  identified  by  the  coordinate  2^  of  its  place  of  symmetry.  It  Is 
sasuoed  that  c(s)  has  a  constants  generally  ccs^lex,  val'%-  c-f  aq  in  each  one  of 
these  subregions.  Tfcj  axial  dimension  of  these  subregions  can  be  taken  smaller 


f 
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in  the  vicinity  of  tbs  active  transducer  to  account  fa-  the  acre  rapid 
decay  of  a(z)  near  the  redieting  surface.  She  z^.- integral  in  the  integral 
equation,  Eq.  H.15,  la  no*  replaced  'ey  a  cusoetiaa  over  the  0  suabregtooc. 

integral  equation  auat  he  eatisfied  at  discrete  points  a  in  the  region  z„: 

r  r 


d  rrfs  -8  i  +  6  —1 

9  '  P  V  pq  a*  f  a 


*«*- 

J  A<va# 


)  da' 


The 


(via) 


VL 


where  6^  is  the  Kronecker  delta  function  which  equals  1  when  p«q  and  zero  when 
p/q.  If  p  is  successively  set  equal  to  ail  valves  of  q,  frot  1  to  Q,  a  set  of  0 


simultaneous  linear  equations  in  cr^  la  obtained. 


In  this  particular  approach 
there  is  as  chrlous  advantage  at  evaluating  separately  the  oyanetric  and  anti* 
sysaetrlc  velocity  components.  The  two  uncoupled  integral  equations,  Sqs.  IV.  12 
now  give  rise  to  two  uncoupled  sets  of  simultaneous  Hn onr  equations.  Hoverer, 
since  the  two  Integral  equations,  Sqs.  IV.  12,  only  extend  over  the  positive  halt 
of  the  region  s_,  the  sane  finite-difference  spacing  results  in  0/2  (instead  ox  0) 
simultaneous  eqt*ations  for  (cc^)s,  and  in  the  saws  nuaber  of  equations  for  fc  )0« 
Thus  conparable  accuracy  la  achieved  by  solving  two  watrlx  equations  e»f  order 
14=0/2,  Instead  of  the  singls  aetrix  equation  of  order  0  which  aust  be  solved  if 
the  velocity  distribution  is  not  split  in  this  fashion.  Tfczse  Matrix  equations 
are  of  the  fora: 


+  2dirg{0i  2d2r|rl’22)  *  ‘  *  £Vg*VZx’ 


1  cul* 


SC*  >“*!> 


2m 


2^r(°) 


[aairi<Vii> 


2*a 


V3i> 

< 

# 

J2. _ 

*-< 

y*s> 

1 

• 

y*. 

',<V 

(VI.2) 


where 


z>L 

F  (s  }  *  7  A  (s  •.,)  4z/ 

YL 


(VI.3) 
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Like  the  Oreen's  function*  of  which  they  a?e  coapoaed,  tbe  function*  r  (t-z') 
here  *  singularity  at  tnf.  tbe  ditgenai  tern*  in  tbe  aatrix  do  not,  nowever, 
display  a  singularity  since  they  are  equivalent  to  *n  integral  of  r(2-a')  over  z' 
which,  like  the  potentials,  1*  a  vsll-btheved  function.  * 

Than  this  aatrix  is  solved  for  tbe  values  of  aQ  tbe  radiation  iapedaoces  SQ 
ere  obtained  in  tbe  font  of  a  gustation 


M  s^+Xi 

fcrja  -  -2(2*a)2top  VYJ  &**/,(*-**)  *  <*asFa(*-*B>3  dt  (w«*) 


VU.  Variational  Solution  of  tbe  Integral  Bauatlan 

As  in  tbe  finite -difference  solution,  it  trill  be  found  advantageous  in  tbe 
present  approach  to  handle  separately  tbe  syamtric  end  entisyasetrlc  components 
of  the  unknown  velocity  distribution.  Our  purpose  In  this  section  is  thus  to 
find  a  rariftticoil  solution  to  the  uncoupled  integral  equations,  Sqe.  17.12,  which 
define  tbe  unknown  velocity  distributions  ssi  or.  fhe  procedure  is  to  con¬ 
struct  functionals  J[a]  which  are  stationary  with  respect  to  first  order  varia¬ 
tions  of  the  soluti<ss  c:g(z)  end  ag(z).  the  steps  deriving  this  functions! 
parallel  those  presented  in  rsf*  J.  Section  7.  Both  sides  of  tbe  Integral  equa¬ 
tion  are  aultiplied  by  a(z)  and  irtegrated  with  respect  to  *  over  the  region  zf+. 
One  thus  obtains  e  relation  between  two  functionals 


B[ag]  -  -  A[agJ 

where 


Sj+L 

“V  ■  / °«WI/ .  ^.(s-z*)  dz#]  4 


(vn.i) 


(VII. 2a) 


3[<*gJ  *  J  j ^(tjir^z-z')  +  6{z-s/)3<^(s/)  dz'  ctz  (VII. 2b 


Sf+  *f+ 


2 

Vasa  we  divide  both  aides  of  8q.  VH.1  by  A  [a]  end  take  their  reciprocal,  wo 


obtain  tbe  following  relation 


ASfa] 


-  A(aJ 


(VIt.3) 


^Subscripts  s  and  a  will  be  csdtted  In  subsequent  .-’uatioss. 
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She  functional  which  displays  tba  desired  Ptationsry  characteristics  is  defined 
by  the  left  side  of  tbs  above  equation 


^  A2fa] 


tm.h) 


lucresseat  &X[a]  associated  vitb  an  incwwat  bet  can  be  vritten  ae$ 


-*  jAJL  _ 

«J[al  -  {B(a}  +  A[os]}  ~i2I  f  Y  }bs{z)  dz'  ds  {VH.5) 

B  !ei  4  vL 

fbr  the  details  of  this  trensfsmstion,  the  reader  ia  referred  to  ref.  2, 

Section  V.  the  functional  J[aJ  is  stationary  if  aj«o.  Since  neither  the  ratio 
2A/B^  nor  the  double  integral  is  identically  zero,  the  cond  jbs  required  to 
sake  the  increaant  ij[s]  vain  «h  thet  tir  *»  Knees  he  -taro.  it  ee»  be 
verified  that  this  aaounts  to  requiring  chat  the  functionals  A[cr]  and  Bfce] 
satisfy  Sq.  VII. 1.  This  equation  is  derived  from  the  integral  equations, 

Eqa.  IV-12.  It  is  therefore  satisfied  by  functions  which  are  solutions  of  the 
integral  equation .  The  functions  which  satisfy  the  integral  equation  therefore 
also  «ka  5J  vanish.  We  hare  tins  shown  that  the  functional  defined  In  *0.  VH.4 
is  stationary  with  respect  to  variations  to  about  the  solution  a(a)  of  the  orig¬ 
inal  integral  equation. 

this  variational  principle  can  be  used  to  obtain  apprcsisRto  values  of  the 
radiation  lapedanee  in  such  the  saee  way  that  the  Bayleigh-Biis  mcitxA  yields  ex¬ 
pressions  for  natural  frequencies  of  Titrating  systems.  The  procedure  veg  prv- 
seated  in  detail  in  ref.  1,  Section  VI  and  will  therefore  only  be  briefly  re¬ 
viewed.  In  the  Rayleigh-Ritz  procedure  a  trial  function  fw  oc(z)  is  constructed 
ia  terns  of  undeteraired  coefficients  associated  with  various  powers  of  *;  It  was 
shown  that  in  the  far  field,  a{z)  be  cf  tbs  fora 


aW  -  for  b  |  Inrgc 

S’  1 


(vn.6) 


In  the  near- fie  Id  of  the  o'.Tay,  for  values  of  s  only  slightly  larger  than  l,  c 
sore  rapid  decay  with  increasing  s  can  be  assessed  to  y isolate  near  field  condi¬ 
tions.  It  is  thus  expedient  to  include  a  tern  proportional  to  s”^.  When  the 
£2pa  between  svrvy  ilasents  are  narrow,  tne  radial  velocity  through  these  gaps 
can  be  taken  equal  to  a  constant,  x  ,  which  of  course  differs  free  gap  to  gap. 
Phase  shifts  between  the  radial  velocities  in  different  gaps  ore  taken  into 
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i 

i 
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scco&t  cy  the  feet  that  the  coefficients  xg  are,  in  general,  complex.  Since 
splitting  the  velocity  distribution  into  synstrie  and  antiayssstric  components 
confines  the  integral  equation  to  the  subregion  sf+,  the  trial  function  need  be 
defined  only  for  positive  values  of  z.  A  trial  function  in  the  following  for* 
is  thus  obtained 


a(z)  =  (-1  +  -§ )  exp(ikz) 

z3  - 


for  z  >  1 


(vn.7») 


a(z)  =•  x  for  z  -b  <  z  <  z  ,  witfe  s >  S>  (VH*7b) 

D  6  o 

(see  Table  is  fi*.  1  for  values  of  *_).  A  refinement  which  will  be  required  when 
tee  gaps  between  the  array  elements  are  large,  involves  taking  into  account  the 
change  of  radial  velocity  within  the  width  of  each  gap.  Instead  of  Mq.  V££.7b, 
we  can  now  assume 


a(s)  »  x  +  x  ,z  for  z  -b  <  z  <  z  +b,  with  z  >  0  .  (VH. 7c) 

go  g±  6  g  6  — 


Experience  with  numerical  calculations,  and  in  particular  ccaqeriscn  with  the 
approaches  described  in  sections  V  and  VI,  will  indicate  whether  the  radiation 
ispedancfc  is  sensitive  to  the  selection  of  the  trial  function  a(z).  Since  the 
resultant  velocity  distribution  is  the  sun  of  two  trial  functions  oo  and  a_,  a 
fairly  nodes t  nusber  of  vr^nom  coefficients  allows  considerable  flexibility. 

As  an  example,  let  us  consider  a  three-elessnt  array  esix  dying  relatively 
esttow  pips.  The  trial  function  in  Eqs.  VII. 7  cow  becomes 

X,  *9 

a(z)  ■*  (-4  +  —)e2p(ihz)  f or  z  >  3L  -r  2b 
zi  z 

a(z)  »  Xj  for  L  <  z  <  L  +  2b  (VU.8) 

By  virtue  of  the  variational  principle  derived  above,  the  best  values  Of  the  un¬ 
known  coefficients  are  those  which  give  «  stationary  value  to  the  functional  J[a]: 

bl[a}/hxn  *  0  with  n  =  1,2,  and  3  (VH.9) 

When  ve  substitute  the  definition  of  J[aj,  Sq.  VH.b,  and  nultiply  all  terns  by 
B(oJ,  these  equations  beecae 

2  A(a]  -  SgStf  jfa]  -  0,  with  n  -  1,1.  and  3  (VH.10) 
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When  ^1*  functionals  A[a]  ere  B[a],  Sqs.  VZI.2  »rs  «v»lnetr4  in  tame  cf 
8q.  YU. 8,  the  former  functional  is  found  to  be  *  linear  function  of the  xsmSSmS 
coefficients  in  a(a),  end  the  letter  e  quadratic  function: 


A[«] 


k 


(vn.ru) 


EI«1  -  ^(bnx *  +  2 

n*l 


•  2Vn  ♦  g  Y*ajk  (vn.llb) 

n 

fbe  coefficients  a  ,b  anl  to _ are  know,  ccaplex  constant*.  When  2qa.  VII. 11 

s  n  SB 

are  substituted  in  the  simultaneous  variational  equations,  &}.  VH.10,  a  set 
of  linear  equations  in  the  undetermined  coefficients  is  obtained: 

3 

(c2  -  1>nJlo|)*n  +  y~* (*-*a  *  *  0  n»l,a,  and  3  (VH.12) 

wfo  " 

These  equations  being  homogeneous,  they  admit  ncn-tririal  solutions  only  if  the 
determinant  of  their  coefficients  vanishes.  One  thus  obtains  a  characteristic 
equation  vhieh  can  be  solved  far  J[a): 


*1  -  bjJla] 

*1*2  “  bl2Jt°J 

V3  *  VtaJ 

Sla2  "  b12J^ 

*2  '  ^»tcJ 

Vi  ■  b23J[aI 

»0 

V3  '  b13J[a] 

V3  '  V[aJ 

*3  *  to^rotj 

(vn.13) 

In  contract  to  the  Raylaigh-Ritz  natural  frequency  calculation,  vhieh  yields  a 
number  of  roots  equal  to  the  order  of  the  aatrix  of  the  characteristic  equa¬ 
tions,  Sq.  VH.13  a daits  only  care  son-vanishing  root,  *  e. ,  cue  value  for  the 
functional  j[s],  fta  reason  is  that  sll  terra  except  those  containing  the  two 
highest  power*  of  j[a]  cancel.  This  result  is  consistent  vlth  the  fact  that  the 
non  -  bosogeneous  integral  equation  admits  only  one  n  on-trivial  .lolutioa.  the 
7alus  of  J(a]  obtsfesd  fro*  S3.  VH.13  is  given  in  ref.  1,  p.  W>,  &j.  VX.12a. 
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Sw  advantage  of  analysing  the  aatieymmtrio  sod  the  symotrie  components  of 
the  Telocity  distribution  separately  la  that  if  M  unknown  coefficient*  are  used  in 
each  of  the  two  trial  functions  thich  describe  the  resultant  Telocity  distribution 
(<y<*a)»  two  deterninante  of  order  K  mat  be  solved.  If  the  Telocity  distribution 
is  not  split  in  this  fashion,  a  determinant  of  carder  at  aust  be  solved  to  achieve 
a  solution  of  comparable  accuracy.  This,  of  course,  is  more  laborious. 

It  is  sham  in  reference  1  that  Then  h/a,  the  transducer  veil  thickness  to 
radius  ratio,  is  sufficiently  email  to  mmfce  the  quantity 

I&-»>  -  (;>'  v  «  1  in  the  oaae  of  amgnetamtrlctiTe  transducers 


|{l-  -  (~)2  ~  «  I  in  the  case  of  piesoelectrlc  transducers  (VU.14) 


the  coefficients  8i  and  8q  can  be  set  equal  to  unity,  fills  in  turn  aekea  th»  ■ 
Green's  function  combinations  T  (z-z4)  ord  «g{z-z~)  in  Xqs.  IV.  11  equal.  With 
this  assumption,  and  making  use  of  the  fact  that  the  functions  r  (x-s*)  are 
symmetric  in  (s-2/)  to  inrert  the  order  of  integration,*  the  functional  in 
Eq.  VH.2a  can  be  written  as 


A“V 


rJ+L  r 

f  l  a  (z')r  (z.z')  az'j  d= 

{  i  1  8 

SJ-L  8f+ 


(vn.15) 


When  q(z)  satisfies  the  Integral  equation,  Eq.  VII. 3  applies,  and  the  functional 
J[a],  2q.  VH.h, equals  -A [a].  If  we  now  compare  Sqs.  VH.15  and  IV. It,  we  see 
that  (Z^,  i.e.,  the  component  of  the  self- Impedance  associated  with  the  radial 
velocity  distribution  (a«a jKt  ),  is  proportional  to  the  sun  of  the  roots,  Jfa  ] 

8  B  0 

and  J(a  J,  of  the  characteristic  equations,  Bq.  VU.13: 

S 

(Z3i)Q  -  top(2sa  f  (J[aaj  4  Jta.3)  (VH.16) 

Thus  for  array  eleaents  whose  ratio  h/a  is  sufficiently  snail,  the  self -impedance 
of  any  eleaent  can  be  obtained  directly  fr css  the  characteristic  equation  without 
previously  having  to  compute  the  undetermined  coefficients  in  the  trial  functions. 

If  the  ratio  h/a  is  too  large  for  this  approaiaation,  the  values  of  the  functionals 


*Tha  expression  "syassetrie”  scans  that  r  and  can  be  interchanged,  i.e.,  that 
r(s-r,,)»r(2,-x).  in  Section  IV,  the  tera  "sysBotric"  refers  to  the  ihet  that 
changing  s  into  •*  does  cot  change  the  value  of  Ts.  P8  i8  aantisymnetric"  in  z, 

'  because  changing  the  sign  of  t  alters  the  sign  of  PQ ,  hut  it  is  sysssetric  in 

j  because  in ter changing  these  two  variables  docs  not  alter  Ta. 

I 
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suet  be  substituted  bade  into  the  corrgaprwxMTig  set  of  M  homogeneous  gqut» 
tloes,  89s.  VH.12,  each  of  which  will  yield  (K-l)  ratios  of  undetermined  coeff¬ 
icients. 

Finally,  the  value  of  one  o?  thss«  urhnova  coefficients,  say  3^,  is  obtained 
by  substituting  those  ratio*  in  the  functional  Ala]  a*  given  by  Sq.  VET. 11a.  As 
noted  earlier,  this  functional  equal*  -J[a],  fer  the  value  of  0(2}  which  satis¬ 
fies  the  integral  equation.  Therefore,  setting  Eq.  VH.lla  equal  to  the  value 
of  J[a]  obtained  from  Tq.  VH.13,  we  can  solve  for  the  remaining  undetermined 
coefficient.  Thus,  for  trial  functions  involving  these  undetermined  coefficient*, 

*i  *  (vn.i?) 

An  alternative  procedure  la  to  select  the  renaming  undetermined  coeffi¬ 
cients  so  as  to  satisfy  the  original  integral  equation  itself.  Unless  the  func¬ 
tional  dependence  of  the  trial  function,  Sq.  VH.8,  on  2  is  the  correct  one,  this 
coefficient  can  not  be  selected  so  as  to  satisfy  the  integral  equation  over  the 
whole  region  zf+.  It  is  advantageous  to  select  the  coefficient  eo  as  to  satisfy 
the  equation  for  a  value  of  2  associated  with  a  relatively  large  value  of  a(z), 
i.e.,  a  value  of  2  which  lie3  close  to  the  active  transducer  element.  In  the 
case  of  a  three-decent  array,  and  for  the  trial  function  assumed  in  Eq.  VU.8, 
one  night  select  x,  so  58  to  satisfy  the  integral  equation  at  z-Ltb.  Previously, 
one  would  substitute  in  the  set  of  homogeneous  equations,  Eqs.  VH.12,  the  value 
of  J[cs]  as  obtained  from  Sq.  VH.13*  One  can  then  solve  for  the  ratios  x^/xy 
The  integral  equation,  Eq.  then  yields 


2,+L 

ri 


-J  a(l*b-z')  dz' 


!i*L 


*3”  5*b  • - - - 3 - - - 

X  3Lt2b  3  U  >  3  {z  ) 


(vn.18) 


Once  again,  tbs  subscripts  "0"  end  "a"  have  been  omitted  in  this  equation 

Having  thus  determin'd  the  coefficients  in  terms  of  which  the  trial  functions 
are  expressed,  the  self-  and  ssutual  impedance  correction  factors,  (zjj)<x  ani 
(Zjk)a,  *re  computed  by  substituting  the  trial  functions  in  Eq.  XV. Ik.  The  far 
field  potentials  can  also  be  obtained,  by  using  the  expressions  given  in  ref.  1, 
Sections  VH  and  VIC. 
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It  nt  sectioned  in  the  introduction  that  the  application  of  the  variational 
principle  to  the  array  is  less  advantageous  thin  its  application  to  the  single 
"aquirfcer."  33se  reason  is  precisely  that,  in  contrast  to  the  self-impedance*  the 
■utual  impedances  cannot  te  obtained  directly  fro*  the  value  of  J{n]  computed  fro* 
Xq.  VU.16,  hat  require  tfcet  Bqs.  VH.12  he  solved  for  the  undeternined  coeffi¬ 
cients  Xg*  It  is  possible  that  a  sore  sophisticated  variational  approach,  whereby 
the  autual  as  well  as  the  self-radiation  iapedanees  can  ho  found,  directly  fro* 
stationary  functionals,  can  he  evolved  hy  formulating  the  array  analysis  ss  a  aul- 
tiple  scattering  problem.  In  this  fonuloti.cn  the  scattering  cross  sections  of 
the  array  elements  would  he  expressed  in  terns  of  stationary  functionals  constructed 
in  accordance  with  the  Levlne-Scfevii^er  vsriatlcnal  forxulstion  of  diffraction 
problene.  The  radiation  iapedunces  would  then  be  related  to  the  scattering  cross 
sections. 
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